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Abstract. We show in spatially one dimensional Madelung fluid that a simple 
requirement on local stability of the maximum of quantum probability density will, 
if combined with the global scale invariance of quantum potential, lead to a class of 
quantum probability densities globally being self-trapped by their own self-generated 
quantum potentials, possessing only a finite-size spatial support. It turns out to belong 
to a class of the most probable wave function given its energy through the maximum 
entropy principle. We proceed to show that there is a limiting case in which the 
quantum probability density becomes the stationary-moving soliton-like solution of 
the Schrodinger equation. 
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1. Madelung fluid: global scale invariant quantum potential 

Let us consider a spatially one dimensional Madelung fluid p] for a single free particle 
with mass m. The state of the system is then determined by a pair of fields in 
position space q as: {p(q), v (q)}, where p(q) is a normalized quantity called as quantum 
probability density and v(q) is velocity field. The temporal evolution of both at time t 
is then assumed to satisfy the following coupled dynamical equation: 
(in 

m— = -d q U, d t p + d q (pv) = 0. (1) 

Here, U(q) is the so-called quantum potential generated by the quantum amplitude 

R = p 1/2 as 

For the case of a spatially one dimensional fluid, the velocity field can always be 
written as the spatial gradient of a scalar function S(q) as 

v(q) = d q S/m. (3) 

One can then use this new quantity to define a complex-valued wave function ip = 
Rexp(iS/h) to show that the Madelung fluid dynamics given in Eq. (JTJ) is equivalent 
to the Schrodinger equation for a single free particle as follows: 

ihd^{ q] t) = -^dl^{q-t). (4) 

Notice that while the equation on the right of Eq. (CQ) is but the conventional 
continuity equation which guarantees the conservation of probability flow, the left 
equation takes the form of Newtonian dynamical equation with a classically absence 
new term appears on the right hand side. In this regards, the term F = —d q U is 
called as quantum force. This fact suggests that the quantum force is the new quantity 
which is responsible for the nonclassical behaviors of Schrodinger equation. Thus, it is 
reasonable to pay serious attention to the the property of the quantum potential. 

Let us mention an important properties of quantum potential that will play 
important roles in our discussion later. First is multiplying the quantum probability 
density with a constant will not change the profile of the quantum potential. Namely, 
the quantum potential is invariant under global rescaling of the quantum probability 
density, namely its own source, as: 

U(cp(q)) = U(p(q)), (5) 

where c is constant. This shows that the quantum potential only cares about the form 
of the quantum probability density and is independent from the strength of the latter 
[2]. It is as if the quantum potential considers the quantum probability density as a 
code in telecommunication system in case of which only the profile of the sequence of 
the binary wave is important, the strength of the received wave itself is of no use. In 
this sense, the quantum potential is of informational nature. For an interesting and 
stimulating discussion concerning this matter see [2]. 
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Among the consequences of the above invariant property is that, first rescaling the 
quantum probability density by the mass of the particle, p = mp will not change the 
dynamical equation on the left of Eq. (JTJ). On the other hand, the continuity equation 
on the right of Eq. ([I]) becomes 

d t p + d q (vp) = 0, (6) 

which can now be read as the equation for the conservation of mass density, rather than 
the conservation of quantum probability density. The other consequence of the scale 
invariance property of quantum potential is that even at points where the strength of 
the quantum probability density is very low, the quantum potential that it generates 
at that point might be very high. In this paper, we shall be interested in a class of 
quantum probability densities with this specific property. 

Next, let us mention another property of quantum potential U(q) that the average 
of the quantum force, F = —d q U, over the quantum probability is vanishing [3] 

'dqd q U(q)p(q) = 0. (7) 



This can be proven easily by assuming p(±oo) = 0. Imposing this into the dynamical 
equation of Eq. (CQ), one reproduces the Ehrenfest theorem [3] 

rnf t = 0, (8) 
where v is the average value of the velocity field defined as v = J dq v(q)p(q). 

2. Local stability and globally self-trapping quantum potential 

Let us now show that local geometrical restriction on the maximum point of the quantum 
probability density, if combined with the scale invariance of the quantum potential will 
determine the global geometrical property of the quantum potential, thus the quantum 
probability density as well. First, since the quantum probability density is vanishing 
at infinity, p(±oo) = 0, non-negative and normalized, it must at least have one local 
maximum point. Let us denote this maximum point by q = Q. It thus satisfies 

d gP \ Q = 0, d 2 qP \ Q <0. (9) 

One first observes that at this point the quantum potential is positive definite 

™~£(-5(¥) ,+ ?)L> a <«* 

Before proceeding, let us write a useful formula for later discussion 



dgP s <9"p 
- s— 



Q p 



which can be shown easily by utilizing the left equation in (jHJ) to be valid for any positive 
integer n. 
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Now, let us put a local restriction on a class of quantum probability densities p(q) 
so that its maximum point stays at the minimum point of the quantum potential U(q) 
which it generates through Eq. (j2J). One therefore imposes 

d q U\ Q = o, d 2 u\ Q > 0. (12) 

Dynamically we are thus looking for a class of quantum probability densities in which 
at least its maximum is temporally stable. Next, let us show that the restrictions given 
by Eqs. (fT2l will uniquely determine the form of U(q) as a function of p(q). First, from 
the global scale invariant property of the quantum potential, then the quantum force is 
also global scaling invariance; so that one has d q U(cp) = d q U(p) for any real constant 
c. It is therefore reasonable to write the quantum force to take the following non-trivial 
form: 

d q U(p) = i(a + ai d q + a 2 d 2 q + a 3 d 3 q + . . .)p s , (13) 

p.s 

where s and dj, % = 0,1,2,..., are arbitrary real number. Evaluating at q = Q and 
using the fact of Eq. (TTTj) one has 

d q U(p(Q)) = a + -( ai d q + a 2 d 2 q + a 3 d 3 + ...)p (14) 

p V 

The left equation in ffl2|) imposes the right hand side of Eq. (}T4"j) to be vanishing. 
Keeping in mind Eqs. ([9]) and the fact that <9™p|q, for n > 3, are fluctuating between 
positive and negative value, d q U\q = can then be accomplished by imposing a = 0, 
cij = for j > 2, and a x is arbitrary, yet non-vanishing. One therefore has 

d q U(p)= ai ^f. (15) 
p s 



Next, let us rewrite Eq. f[T5"|) as follows 



d q p _ d q p 



d q U (p) = aid q In p s = aisd q lnp = ais—^— = a s -^—, (16) 

where we have denoted a s = a\s. Now, taking spatial derivation on both sides of the 
above equation and using the left equation in Eq. ([9]), one gets 



dp 



d qP /<V\2 8 2 p 

a s ^— — a.A ) _ = a.,^— 

p Q 



V p J Q p 

Comparing this fact to the right inequality in ( fl~2l) and keeping in mind the fact that 
2 p\q < 0, one concludes that a s must be non-positive. One can then verify that 
any quantum probability density that satisfies Eq. ( TT6l) satisfies all the requirements 
that we set at the beginning. Moreover, assuming p(±oo) = 0, Ehrenfest theorem is 
automatically satisfied 

J dq pd q U = a s J dq d q p = 0. (18) 

To proceed, for simplicity of notation, let us rewrite Eq. (TIB"]) as follows 

dM = -a^-, a>0. (19) 
P 
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It can be readily integrated to obtain 

p(<?; a ) — ~77r\ ex P ( — ~U (q; a)] , (20) 
Z{a) v a * 

where Z(a) = J dq exp(— U/a) is a normalization constant independent of q. We 
shall show later that U(q) can be interpreted as internal energy density. Bearing 
this in mind, then the quantum probability density given in Eq. (TSUI) resembles in 
form with the Maxwell-Boltzmann-Gibbs (MBG) canonical distribution in equilibrium 
thermodynamics. It is thus suggestive to apply thermodynamics formalism to further 
study the property of quantum probability density given in Eq. (1201 [5]. 

Next, let us recall that in quantum mechanics p(q) gives the essential information 
on the position of the particle [21 El [7]. In the so-called pragmatical approach of quantum 
mechanics, p(q) is given meaning as the probability density that the particle will be found 
at q if a measurement is performed. On the other hand, in the ontological approach, 
p(q) is argued as the probability density that the particle is at q regardless of any 
measurement. It is thus reasonable to quantify the randomness encoded in p(q). One 
obvious way is then to use the differential entropy or the so-called Shannon information 
entropy [8] over the quantum probability density given by 

H[p] = -J dqp(q)\np(q). (21) 

It gives the degree of localization of the wave function in position space. 

One can then show that the canonical quantum probability density of the form 
given in Eq. (T2Q|) maximizes the Shannon entropy provided that the average quantum 
potential is given by [9] 

U = J dqU(q)p(q). (22) 

Hence, the quantum probability density developed in the previous section satisfies the 
so-called maximum entropy principle [10] . It has been argued that the maximum entropy 
principle is the only method to infer from an incomplete information, which does not 
lead to logical inconsistency [TT]. The self-trapped quantum probability density can then 
be seen as the most probable quantum probability density given its average quantum 
and kinetic energy [12] . 

Combined with the definition of quantum potential given in Eq. ([2]), Eq. (!20j) 
comprises a differential equation for p(q) or U(q) subjected to the condition that p(q) 
must be normalized, / dqp(q) = 1. In term of quantum potential, one has the following 
nonlinear differential equation 

%v-±{d 9 u)*-^u = o. (23) 

Figure [1] shows the solution of Eq. ff23l) with the boundary conditions: d q U(0) = 
and U(0) = 1 for a = 1. All numerical solutions in this paper are obtained by putting 
m — h — 1. The quantum potential is shifted down so that its global minimum is 
vanishing. One can first see that the maximum point of the quantum probability density 
and the minimum point of the corresponding quantum potential coincide, thus satisfies 
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our requirement. Yet, what makes even interesting is that, though we only requires 
the quantum potential to trap the area of the quantum probability density around its 
maximum, it turns out that the resulting quantum potential is convex everywhere and 
becomes the global trapping potential for its own source: quantum probability density. 




Figure 1. The profile of quantum probability density and its corresponding quantum 
potential which satisfies Eq. (f!?5]) . 

Next, one can also see that the solution plotted in Fig. [1] possesses blowing-up 
points at q — ±q m , namely U(±q m ) = oo [I]. Let us first prove that the blowing-up will 
certainly occur at finite point from the origin as along as U(0) = X is not vanishing. 
To do this, Let us define a new variable u = d q U. The nonlinear differential equation 
of Eq. (123]) then transforms into 

_ 1 2 4ma T7 , nAS 
d qU = —u 2 + —^U, (24) 




Figure 2. u(q) and u(q). See text for detail. 
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Moreover, the boundary condition translates into u(0) = d q U(0) = 0. Let us now 
consider the following nonlinear differential equation 

1 _ 4ma , , 

^ = Ta" + ~W < 25) 

where X = U(0) with u(0) = 0. Since U(q) > U(0) = X, then it is obvious that 
\u(q)\ > \u(q)\. 

On the other hand, one can solve the latter differential equation of Eq. (1251) 
analytically to have 

u(q) = dtan(gq), d = — V2m, g = —\ / 2mX. (26) 

It is then clear that at q — ±q m = ±n/(2g), u is blowing-up, namely u(±q m ) = ±oo. 
Recalling the fact that \u(q)\ > \u(q)\, then u(q) is also blowing-up at point q = ±q m , 
u(±q m ) = ±oo, where q m < q m . See Fig. [2j Putting this into the original nonlinear 
differential equation of Eq. (1231) . one concludes that U (q) is also blowing up at q = ±q m , 
U(±q m ) = oo. Notice that even though u(—q m ) = —oo is blowing-up to minus infinity, 
U(—q m ) = oo is obviously blowing up into positive infinity. Next, it is clear that the 
blowing-up is due to the existence of the nonlinear term on the right hand side of 
Eq. ( l23l) . Hence, finally one can safely say that for any non- vanishing U(0) = X, 
the corresponding quantum probability density possesses only a finite range of support, 
q £ [— <?m,<?m]- The case when U(0) = will give the trivial solution U(q) = for the 
whole space q so that p(q) is unnormalizable. The above fact also confirms our assertion 
in Section I that the quantum potential might take large value even at points where the 
corresponding quantum probability density is very small. 

In Fig. [3] we plot the blowing-up point q = q m , namely half length of the spatial 
support of the quantum probability density against the value of the quantum potential 
at the global minimum: U(0) = X. One observes that q m is decreasing as we increase 
X for fixed a = 1. This can be understood directly from Eq. ( 1261) . One can also confirm 
that the occurrence of blowing-up is the case only when 1^0, namely linix^o Qm — °°- 



X 



Figure 3. The half length of support q m plotted against the variation of the global 
minimum t/(0) = X . 
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3. Solitonic wave function 

Now let us proceed to study the behavior of the quantum potential as one varies the 
non- negative parameter a. Figure 0] gives the variation of the blowing-up point, q m , 
thus the size of the range of the support against the variation of the parameter a. 
This is obtained by solving the differential equation of Eq. (1231) with fixed boundary 
conditions: d q U(0) = and U(0) = 1. One first observes that as a is increased, 
q m decreases and eventually vanishing for infinite value of a. This shows that the 
quantum probability density is becoming narrower for larger a while kept normalized; 
and eventually collapsing onto Dirac delta function for infinite value of a. 




Figure 4. The half length of the support against the variation of a. 

A very interesting phenomena is observed as one decreases the parameter a toward 
zero. One finds that the blowing up point q m is increasing and eventually converging 
toward a finite value go for a = 0, 

limg m (a)=g - (27) 

This suggests to us that the quantum potential and the corresponding quantum 
probability density are also converging toward certain functions for vanishing value 
of a: 

\imU(q;a) = U (q), lim p(q; a) = p (q). (28) 

a— >0 a— +0 

Let us discuss this situation in more detail. In Figure [5] we plot the profile of the 
quantum probability density and the corresponding quantum potential for several small 
values of parameter a with fixed boundary conditions: d q U(0) = and U(0) = 1. One 
can then see that as a is decreased, the quantum potential is becoming flatterer inside 
the support before blowing-up at q = ±q m (a). One might then guess that at the limit 
a = 0, the quantum potential is perfectly flat inside the support and is infinite at the 
blowing-up points, q = ±qo- Let us show that this guess is correct. To do this, let us 
denote the assumed constant value of the quantum potential inside the support as U c . 
Recalling the definition of quantum potential given in Eq. (T5]), one has 

- 7^o(<7) = U c Ro(q), (29) 
2m 



1 II 

where Rq = p . This has to be subjected to the condition that i?o(±9o) = 0. 
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a = 0.5 



a = 1 



oo 




Figure 5. The profile of self-trapped quantum probability density and its 
corresponding quantum potential for several small values of parameter a. We also 
plot the case when a = obtained analytically in Eq. (|30|) . 



Next, solving equation (T2"9~j) one obtains 
R (q) = A cos(k q), 



(30) 



where A is a normalization constant and ko is related to the quantum potential as 



The boundary condition imposes k q Q — tt/2. In Fig. [5j we plot the above obtained 
quantum probability density, po(q). One can see that as a is decreasing toward zero, 
p(q; a) obtained by solving the differential equation of fl23l is indeed converging toward 
Po(<?) given in equation fl30|) . This confirms our guess that at a = the quantum 
potential takes a form of flat box with infinite wall at q = ±go- 

Let us now take {po(q) , v (q)} as the initial state of the dynamics. Here v (q; 0) = v c 
is a uniform velocity field with non- vanishing constant value only inside the support. 
Since at a = the quantum potential is flat, then inside the support the quantum force 
is vanishing: F = —d q U = 0. Inserting this into the dynamical equation of Eqs. (CQ), one 
has dv/dt = 0. Hence the velocity field at infinitesimal lapse of time, t = At, remains 
constant and uniform. This in turn will not change the initial probability density, but 
shift it in space by Aq = v c At: p(q; At) = p (q — v c At; 0). Accordingly, the support is 
also shifted by the same amount. This will repeat in the next infinitesimal time lapse 
and so on and so forth so that at finite lapse of time t, both the initial velocity field and 
quantum probability density remains unchanged but are shifted by an amount Aq = v c t. 
One thus has 



where we have put uj = k v c . 

Before proceeding, let us give physical meaning to the average quantum potential. 
To do this, let us calculate the ordinary quantum mechanical energy given by (E) = 




(31) 



P(Q] t) = Po(q - Vet] 0) = Al cos 2 (k q - u t), 



(32) 
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Jl° qo dq '4 ) *(q)( — h 2 /2m)d 2 ip(q) . Writing the wave function in polar form one has 

(E) = r dq ( - ^-Rd 2 q R + ^-R 2 (d g S) 2 
J-qo v 2m H 2m 

--Rd q Rd g S - —R 2 d 2 S). (33) 
m q q 2m q ' y 1 

The first term on the right hand side is equal to the average quantum potential, 
U = J dqUp. Next, defining kinetic energy density as K(q) = (m/2)v 2 (q), the second 
term is equal to the kinetic energy K = J dqKp of the Madelung fluid, which for our 
stationary state is given by K = (m/2)v 2 . Further, for a uniform velocity field, the last 
term is vanishing, d 2 S = md q v c = 0. Again for a uniform velocity field, since R(q) is 
an even function and d q R(q) is an odd function then the third term is also vanishing. 
Hence, in total, the quantum mechanical energy of a self-trapped wave function for 
a = moving with a uniform velocity field can be decomposed as 

(E) = U + K. (34) 

One can then conclude that Uq must essentially be interpreted as the rest energy of the 
single particle. Namely it is the energy of the particle when it is not moving so that 
K = 0. Moreover, since inside the support the quantum potential is flat given by U c , 
one has Uq = fl° Uo(q)p(q) = U c . Recalling Eq. (T5T|) . one finally obtains 

{E) = If + \ mv '' (35) 

Let us now give the corresponding complex- valued stationary wave function ip(q). 
One thus needs to calculate the quantum phase S which can be obtained by integrating 
d q S = mv c to give: S(q;t) = mv c q + cr(t), where a(t) depends only on time. One 
therefore has ip s t(q',t) = A cos (k (q — v c t)) exp (^(i/h)(mv c q + cr(t))j, where q G M. t = 
[v c t — qo, v c t + q ]. Inserting this into the Schrodinger equation of Eq. (jlj) and using Eq. 
( 135]) one has a(t) = —(E)t modulo to some constant. Putting all these back, one finally 
obtains the following solution: 

ip st {q;t) = A cos {k (q-v c t)) exp ^(mv c q- (E)tfj, (36) 

where q G M. t - 

Equation ( |36|) is of soliton type. This suggests a direct association of the wave 
function to a particle by considering the wave function as a physical field. In this regard, 
the continuity equation of mass density of Eq. ([6]) is becoming relevant. On the other 
hand, we have also shown in the previous section that the localized-stationary-traveling 
solution belongs to a class of wave function which maximizes Shannon entropy, which 
suggests that it is a probabilistic wave field. Hence, one arrives at one of the old problem 
of quantum mechanics concerning the physical status of the wave function. 

4. Conclusion and Interpretation 

By exploiting the scaling invariant property of the quantum potential, we show that 
the requirement of local stability on the maximum of quantum probability density leads 



11 



us to a class of quantum probability densities which is globally trapped by its own 
quantum potential with finite-size spatial support. It turns out that they belong to a 
class of wave function which maximizes Shannon entropy given the average quantum 
potential. Further, we show that for a single free particle quantum system, there is 
an asymptotic limit in which the self-trapped wave function is traveling while keeping 
its form unchanged. This fact thus suggests to us to associate the localized-traveling 
quantum probability density as a real particle. In contrast to this, in conventional 
formalism of quantum mechanics, one usually choose a plane wave to represent a free 
single particle. 

In our formalism of a single particle as a localized and self-trapped wave function, 
we showed that a particle possesses an internal energy which is absence if one use a plane 
wave instead. It is equal to the quantum potential. On the other hand, we showed in 
the beginning of the paper that the quantum potential is invariant under the global 
rescaling of its own source, namely the quantum probability density. It depends only 
on the form of the latter. It is a surprising fact, since, usually energy is an extensive 
quantity with respect to its source. Hence, the internal energy is of different nature 
from the ordinary one. Since it only records the profile of its own source, one might 
conclude that its nature is informational, rather than material. 

An interesting point is left unexplored. On one hand, we showed that the canonical 
form of quantum probability density given in Eq. ( 1201) is the consequence of the scale 
invariant property of the quantum potential. On the other hand, we also showed that 
it maximizes Shannon information entropy given its average quantum potential. One 
may then expect that these two facts are related in a nontrivial way. 
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